Current induced distortion of a magnetic domain wall 
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We consider the spin torque induced by a current flowing ballistically through a magnetic domain 
wall. In addition to a global pressure in the direction of the electronic flow, the torque has an internal 
structure of comparable magnitude due to the precession of the electrons' spins at the "Larmor" 
frequency. As a result, the profile of the domain wall is expected to get distorted by the current 
and acquires a periodic sur-structure. 
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With the advent of " spintronics" , which aims at using 
the spin of charge carriers in devices, electronic trans- 
port in ferromagnets is being revisited from a different 
viewpoint. The focus has been for a long time on the ef- 
fect of magnetism on transport properties (e.g. Magneto 
Resistance, MR) but it is now realized that the elec- 
tronic current can be a tool to change the magnetization 
direction. The relevant effect, known as spin torque has 
attracted considerable interest recently in the context of 
ferromagnetic-normal metal-ferromagnetic trilayers 0- 
There, the first magnetic layer acts as a spin filter, and 
the incident polarized electrons exert a torque on the sec- 
ond layer. At the heart of this physical effect is the fact 
that spin currents are not preserved when electrons cross 
a magnetic layer, and as a result some angular momen- 
tum is transferred to the magnetization. It was shown 
experimentally that for strong enough current densities 
this mechanism can lead to magnetic reversal Q. This 
demonstrates the feasability of current controlled mag- 
netic memory cells, but the current needed for complete 
reversal might be too high for industrial implementation. 
An alternative would be to use the current to move a do- 
main wall (DW) in between two stable positions. The 
idea that a current can apply a force on a DW is due to 
Berger 4] in the seventies. It is the aim of this letter to 
study in detail the spin torque exerted on a DW in the 
presence of an electric current. Our main finding is that 
in addition to a global pressure, the torque has a spatially 
dependant component that will lead to a deformation of 
the DW in a periodic sub-structure. 

At the root of understanding the spin torque in a DW is 
the question of what happens to the spin of a conducting 
electron when going through a DW. Two extreme cases 
can be considered called " interface" and " adiabatic" . A 
very sharp domain wall (expected in constrictions for ex- 
ample) , can be treated as an interface on which the elec- 
trons can be specularly reflected f^. This leads to the 
giant magneto resistance effect as a magnetic field will 
remove the DW and the associated extra resistance. Be- 
cause the electronic spin is conserved during this process, 
no spin torque is exerted on the DW. On the other hand, 
in a very long domain wall, the electron's spin will adi- 
abatically adapt itself to remain aligned with the local 



magnetization under the effect of the "Larmor" preces- 
sion 1^. There, no MR is observedj^Q but each electron 
going through the DW will flip its spin and give a quan- 
tum h of angular momentum to the wall inducing a global 
pressure. 

Experimentally, the wall resistance (DWR) has been 
measured in macroscopic systems where it is of the order 
of a few percent per wall. This resistance results from 
a slight mistracking of the conduction electrons' spins 
which mixes the majority and minority channels within 
the wall 0, H, ■ Numerous measurements have been re- 
ported in the last few years finding both positive |^ Q| 
or negative [10| effects. The spread in the experimen- 
tal results probably reflects how difficult it is to extract 
the DWR among other contributions in series coming 
from domains (like the Anisotropic MR) . Other (ballistic) 
models were developed to explain negative effects, 

but it was later recognized that proper band structure 
calculations are needed to get reliable quantitative 
results. Nevertheless, recent clean experiments 0| in 
nanostructures of perpendicularly magnetized materials 
(FePd) have demonstrated that the spin scattering mod- 
els developed in |^ can account for the MR in domain 
walls in 3d metals. Concerning torque effects, very re- 
cently several groups have been able to push DWs with 
a current [isj . 

In this letter, we proceed as follows. First, an heuris- 
tic argument is given to explain the origin of the spatial 
structure of the torque. Then, we introduce our (ballis- 
tic) model and point out that some key features of the 
band structure must be taken into account. Finally we 
calculate the spin torque and evaluate the corresponding 
distortion of the DW proflle. 

Heuristic argument: Let us follow an electron going 
through the (Neel) DW of size A^, sketched in FigQ] (a 
Bloch wall is completely equivalent for MR and torque) . 
The process is described in Fig|5] Before the wall, the 
electron's spin is aligned with the local magnetization m 
(a). Then, when entering the DW, m begins to rotate, 
and a small angle a starts to build up between the spin 
and the local magnetization (b). As soon as a is not zero, 
the electron's spin starts to precess around the direction 
of an effective magnetization with a period equal to the 
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FIG. 1: Schematic of the domain waU of length An, and 
transverse dimension d. The arrows stand for the magne- 
tization direction m that makes an angle 8{x) with the z- 
axis. The same x, y and z axis are used for both the 
space and spin basis for convenience. The rotating frame 
(u, V, w) — (drin/dO, drh/dO x in, m) follows the magnetization. 



"Larmor" precession length Xl ((c) and (d)). At the end 
of every period, the electron's spin is back onto the local 
magnetization direction At x — Al/2, a reaches a 
maximum. Over this distance, m has rotated by an angle 
ttXl/2Xw which thus gives an upper bound for a. Hence, 
V = Xl/Xw is the parameter that controls the crossover 
from the adiabatic limit (77 <C 1) to the interface one 
('7»1)- 
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low, it gives rise to the spatially dependant part of the 
torque inducing a deformation with a period A^. 

Model. Let us introduce a simple model for the con- 
ducting electrons through the DW of FigQ] The system 
is a ferromagnetic wire running along x and of typical 
transverse dimension d. We approximate the wall by a 
linear rotation of the local magnetization on a scale A^. 
The unit vector 771(2;) lies inside the xz-plane and makes 



an angle 0{x) = nx/Xw 
The Hamiltonian reads, 

H = -- 



with the z axis inside the DW. 



2771* 



— m[x) 



(1) 



Here Jcxc is the exchange energy, m* the effective mass, 
= (fx, CTj,, (Tz) the vector of pauli matrices and the 
Fermi energy is noted Ep. The spin dependant part of 
H reads, 



777(0;) • a = 



x <Q 



9{x)i 

Atm <C X 



< X < Xy 



(2) 



with the rotation matrix defined as Rg^x) — e^"^''^^^-'/^. 
In the region outside the wall, the eigenstates of H are 
plane waves in both the transverse j/z-plane, with to- 
tal momentum k^, and in the a;-direction, with mo- 
mentum fc| (fc|) respectively for majority and minor- 
ity electrons. In the transverse direction, is quan- 
tized in units of 27r/d and for the corresponding eigen- 
statc to be propagating in the x-dircction, one needs 
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yf2m*/h?{EF ± Jcxc/2) - (k^Y > which leads 
to iV| {Ni) propagating channels for the majority (minor- 
ity) spin. Once these propagating channels are defined, 
the natural way to calculate physical quantities such as 
the current or the spin current is to use the Landauer- 
Buttiker theory and fill up the different eigenstates of the 
system. This has been widely used in the literature and 
eventually leads to the Landauer formula for the conduc- 
tance g = e^/hTi ti^ where t is the transmission matrix 
of the system. 







FIG. 2: Cartoon of the effect of the Larmor precession al- 
lowing the electron spin to follow the local magnetization m. 
The thin arrow represents rh while the thick one stands for 
the conducting electron spin. 

The (almost) adiabatic case is the most interesting re- 
garding spin torque because the incident spin up electron 
ends up with the (almost) opposite spin after crossing the 
wall, losing an angular momentum h in the process. In 
return, since the total angular momentum is conserved, 
a torque of equal magnitude and opposite sign is exerted 
on the magnetization within the wall. This is the global 
pressure mentioned earlier. The spin precession around 
the slowly rotating effective magnetization does not con- 
serve the angular momentum either, and, as detailed be- 



E 

(a) 



(b) A 



Maj min 



Maj min 



FIG. 3: Cartoon of the position of the majority and minority 
channels as a function of the perpendicular energy. The num- 
ber of majority channels is N] — Na + Nt- (a) simple picture 
where Nb = Ni. (b) "realistic" magnet where Nb ^ Ni. 

In this simple picture, the iVj^ down channels are ex- 
actly matched by the first = up channels with 
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equal perpendicular energy, and the polarization of the 
current is entirely due to the remaining Na = -/V| — iVt, = 

~ channels that have large perpendicular en- 

ergy (see Figl^l (a)). Hence their longitudinal energy 

\J Jexc/rn*), is low and the corresponding minority 
channels may not propagate. This stems from the as- 
sumption that the up and down bands have exactly the 
same shape, which is obviously erroneous. It was recently 
clearly underlined by Mazin |lq that by oversimplifying 
the band structure of magnetic metals, one introduces 
symmetries in the system that ought not to be there. In- 
deed, when the polarization of the current Pi is taken 
equal to the polarization Pn i.e.. 



Pj = 



h-h 



= Pn^ 



Na 



(3) 



this artificial symmetry can lead to erroneous physics. 
For instance, after an "interface" wall, all the channels 
with high perpendicular energy are blocked (the Na = 
(7V| — iV|) in Fig|31(a)) and the polarization of the current 
drops to almost zero. In a real system however (sketched 
in Fig. ini(b)), the number Na of blocked channels would 
not be equal to (iV| — iV^) and the system would retain 
some polarization P/ « {Nii — Ni)/{N'^—Ni). A complete 
treatment of the problem would require ab-initio calcu- 
lations as described in 16]. Here, we break this artificial 
symmetry by giving a different weight to the different 
channels in the Landauer formula, i.e. allowing Pj and 
P/v to be independent numbers. This can be viewed as 
replacing the ferromagnet by a non-magnetic reservoir in 
series with a normal-ferromagnetic interface that acts as 
a spin filter. It leads to the following modified Landauer 
formula, 



Pj-P> 



N 



1-Pn 



(4) 



where the CTz matrix applies in spin space and only to 
those channels that can propagate for both majority and 
minority spins. We emphasize that we introduce this 
modification to take into account the band structure of 
the magnet and do not question the Landauer-Buttiker 
formalism itself. This correction is crucial here since, as 
will be seen later, by symmetry the spin torque vanishes 
when P/v — Pj. 

For a given spinor wave function ^'(x), t he sp in current 
flowing along the x-direction is defined as 
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dx 



(5) 



This is not a conserved quantity, and the corresponding 
loss of spin current is identified with the torque y/fl] . The 
global torque is then Ttot = Js(a; = 0) — Js{x = Xw) while 
the local torque (per unit distance) is. 



d_ 

dx 



Js{x) 



(6) 



From the last equation, we derive for the torque per unit 
voltage:. 



dV 



e_d_ 

Att dx 



Re Tr 



1 



Pi -Pi 



N 



l-P 



N 



t{x)aV{x) (7) 



where the generalized transmission matrix t{x) gives the 
amplitudes of the different modes inside the wall. 

Domain wall close to the adiabatic limit. We now pro- 
ceed with the treatment of Eq. ^ and Eq. Q in the limit 
of a long wall. The first step consists in writing the 
schrodinger equation in the basis aligned with the local 
magnetization. An eigenstate "^{x) is then written as. 



^-(2;) = Reia^)<S>ix) 



(8) 



An effective equation is obtained for ^{x) (see Q fo^' ^ 
similar treatment, the exact solution of the model can 
be found using the "spin spiral" state, see The 
solution is then formally expanded in series of (l/A^,) and 
we match ^'(a;) and its derivative at x — and x — Xw 

In first order in (l/QkXw) where Qk ~ 
for incoming majority electrons in the wall: 



fc|fc|, we get 



$l(x) = ^-p= and ^i{x) 



'iQkXy 



(9) 



Pfc 
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PJ-e'^ = 



-ik' 



and similar solutions for incoming minority spins. The 
wave functions are normalized to carry unit fluxes and 
only the longitudinal part has been written. We have 

defined Pfe = (A:}-A:;)/(fcjl+fcJ), and^ = 6'^''"+''"^^™. We 
point out that the oscillatory first term in the brackets is 
the expression of the Larmor precession while the second 
one is a reflected wave. In the region x < the reflected 
wave takes the form. 



4 QkX, 



(10) 



Spin torque in a macroscopic DW. To proceed with the 
calculation of the conductance and torque, one incorpo- 
rates the expression for the wavefunctions into Eq. and 
Eq. ^ and performs the average over transverse momen- 
tum. This average is done in two steps, first over ^, which 
is taken to be a random phase and then over k±. Also, 
the Fermi wave length A p and Larmor precession length 
Xl are defined as. 



Ap — 2tt\ 



2m* Ep" 



3l 



(11) 



We get for the correction to the conductivity due to 
the presence of the DW,Ag/g = -Pa,A|/(64A^). This 
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correction is (too) small because in a ballistic model, the 
calculated quantity is the reflected part of the wavefunc- 
tion due to the potential step. Diffusive models, on the 
other hand, neglect this contribution and estimate the 
resistance due to spin mixing between the up and down 
electrons. It turns out that in the macroscopic case, the 
latter dominates. In the interface limit however, the bal- 
listic contribution becomes important. We calculate the 
torque per current along the local frame {u,v,w) from 
the wavefunctions expressions as. 



dl 



PTy (X) 

dl 



e 2Xu 



Pi 



{Pi - Pn) cos [2n— 



- — (Pj-P^) Sin (2TT— 



(12) 



Equation (|12|l is the central result of this letter. The 
torque consists of two terms. The first one, proportional 
to Pj pushes the wall in the direction of the electrons 
and does not depend on x. The second part of Ea. lfT^ is 
much more interesting since it leads to a deformation of 
the wall on the scale of A^. Note that although its net 
contribution to the global torque is small, its intensity is 
of the same magnitude as the first term. However, the Na 
channels where the perpendicular energy is high (hence 
the down electrons cannot propagate) do not contribute 
to the torque since no precession is possible inside those 
channels. As a result, this component of the torque van- 
ishes when Pj = P]y, hence the necessity of taking into 
account the fact that P/ 7^ Pn in real materials. We 
point out that although our calculation is done for a bal- 
listic model, the result should hold for realistic domain 
walls. Indeed, both the spin diffusion length (about 50 
nm for Ni) and the mean free path (a few nm) are larger 
than Al (of the order of 3 nm for Ni) and at the Lar- 
mor length scale, the electrons can thus be considered 
ballistic. In a system where the mean free path would 
be smaller than A^ but still with a (rather) large spin 
diffusion length, our conclusions would remain qualita- 
tively correct, though a quantitative treatment might be 
needed. 

In order to quantify the expected deformation of the 
wall, one would have to solve self-consistently the torque 
equation together with the reaction from 771(2;) linked to 
the "stiffness" of the wall (using for instance the Landau- 
Lifshitz- Gilbert equation). Here we simply evaluate the 
magnitude of the deformation by comparing the energy 
given to the wall via the torque to the total wall energy. 
Typical quantities for Ni are A^ = lOOnm, Xf = 0.2nm 



, Xl = 3nm, and Pj — Pn ~ 0.5. A current density 
j = 10^°A/m^ would then leave in a Ni wall 10"^ J/m~^, 
which is one tenth of the wall surface energy. Hence the 
distortion in angular gradient should be of the order of 
several %. This is a significant effect considering that the 
chosen current density is one order of magnitude below 
that used in tri-layer spin torque experiments. 

Constriction "interface " limit. In a constriction (small 
transverse direction d), the DW width is expected to 
scale with d [isL |20| and the system could be driven to 
the interface limit. As explained earlier, the spin of the 
transferred electrons would remain mostly un-flipped and 
the torque would decrease as a result. Hence, while in 
this regime the MR gets larger, the torque described in 
this paper would drop. In addition, since in this limit 
the reflection is higher, another contribution to the to- 
tal pressure could come from the transfer of momentum 
from the reflected electrons to the DW. Although it is 
not clear to us whether this change of momentum would 
actually push on the wall or on the atoms, its effect would 
be much smaller than that in a long wall. 

Conclusion. The torque generated by a current on a 
domain wall in a ferromagnetic metal has been studied 
in unconstrained DWs. It is composed of two contribu- 
tions of comparable magnitude. The first one is a global 
pressure resulting from the loss of angular momentum 
of electrons crossing the wall. The second is due to the 
precession of electrons' spins inside the wall which gen- 
erate a periodic torque. The first effect can be used to 
move the DW with a current. Depending on pinning, 
the current density necessary to dislodge the DW could 
be smaller than that necessary to reverse a full magne- 
tized layer in a spin valve. The effect can then poten- 
tially be useful in spin electronic devices where a DW 
switched between two stable positions could be used in 
the gate of a three terminal device. Moreover, the pe- 
riodic torque will distort the wall's internal structure in 
a significant manner when current densities of the order 
of lO^^A/m^ are driven through it. This may help the 
depinning process and might also be able to switch the 
wall between different types (Bloch and Neel for example) 
in a similar manner to the predicted temperature effect 
in a constriction The distortion may be measured 

by several techniques including polarized neutron reflec- 
tivity (PNR) and domain wall resonance. We also infer 
that a current flowing parallel to the wall would produce 
a surstructure of the same period, which may be easier 
to measure with PNR. 
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